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Ground states in relatively bounded quantum 
perturbations of classical lattice systems 

D.A.Yarotsky * f 



Abstract. We consider ground states in relatively bounded quantum 
perturbations of classical lattice models. We prove general results about such 
Q \ perturbations (existence of the spectral gap, exponential decay of truncated 

correlations, analyticity of the ground state), and also prove that in particular 
the AKLT model belongs to this class if viewed at large enough scale. This 
immediately implies a general perturbation theory about this model. 

Key words: ground state, relative boundedness, AKLT model, cluster 
Tj- ' expansion. 
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^ ■ 1 Introduction and results 

o ' 

^ ! It is generally expected that if a ground state of a quantum lattice system 

^| is in a non-critical regime characterized by the presence of a spectral gap 

J ; and exponential decay of truncated correlations, then the system remains in 

this phase under sufficiently weak perturbations of a general form. Relevant 
rigorous results are now available in the case of weak quantum perturbations 
of some classical models |U IHl CE1 HZl UHl E21 ESJ EHl ED] - Most of these results 
concern perturbations which are bounded and small in the norm sense. How- 
ever, Kennedy and Tasaki obtained in general results for perturbations, 
which are only relatively bounded, in some special sense, w.r.t. the classical 
Hamiltonian. Moreover, using a special transformation of the Hamiltonian, 
they applied this perturbation theory to the dimerized AKLT model, which is 
a genuinely quantum S , f/(2)-invariant model. The type of relative bounded- 
ness they used does not, however, seem to allow an extension of their result to 
the non-dimerized, fully translation invariant case. In this paper we consider 
perturbations relatively bounded in the quadratic form sense, which appears 
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to fit naturally in this and some other contexts. We prove general results for 
gapped classical models with a simple ground state and then apply them to 
the non-dimerized AKLT model. 

We consider a quantum "spin" system on the lattice Z". Throughout the 
paper we consider only translation invariant interactions. Each site x € Z" 
is equipped with a Hilbert space 7i x , possibly infinite-dimensional. In the 
sequel we use the notation 

for Hilbert spaces, corresponding to finite subsets of the lattice. We assume 
that Ti x has a preferred vector denoted fl x . The corresponding product state 
will be denoted by ^a,o : 

Also, we fix some finite set Ao C Z", which will be the interaction range. The 
(formal) Hamiltonian has the form 

H = H + $, 

where H is the classical part and <3> the perturbation. The classical Hamil- 
tonian Hq is given as 

Ho = ^ h x . 

xez" 

Here h x is a self-adjoint, possibly unbounded operator acting on H.a +x, where 
Ao + x is a shift of Ao- The Hamiltonian Hq is classical in the following sense. 
If 7i x is finite dimensional, then we assume that in each 7i x there is an 
orthogonal basis containing Q x and such that the product basis in 7Ya ()+x 
diagonalizes h x . We extend in a natural way this assumption to the case of 
infinite dimensional 7i x by assuming that for each H. x an orthogonal partition 
of unity, containing the projection onto fl x , is given, and h x is a function of 
the product partition in Ha 0+x . Furthermore, we assume that Vl\ + X is a 
non-degenerate gapped ground state of h x : 

K^Ao+xfi = 0, h x \ HAo+x QQ Ao+x0 > 1. (1) 

Now we describe the perturbation. It is given by 

$ = E **> 

x& v 

where 4> x is a (possibly unbounded) symmetric quadratic form on H\ 0+x , 
bounded relative to the quadratic form corresponding to h x , i.e. the domain 
of <p x contains Dom(/z^ 2 ) and 

\(j> x {v,v)\ < a\\h x /2 v\\ 2 + (3\\v\\\ j.GDom(ftf) (2) 
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with some a, (3. We assume that a < 1. The form (j) x actually need not be 
closed and generated by an operator, though in all examples we consider it 
is. 

If A is a finite volume and 

#a,o = ^2 hx > $A = X] 

3i:Ao+iCA x:Ao+xCA 

then $a is again bounded relative to Ha,o with the same a, because, clearly, 
Dom {H^Iq) C Dom ($a) and, by adding up (J2J), 

|$aMI < a||^ \;f + |A|/?|M| 2 , v G Dom(^ 2 ). 

It follows from the KLMN theorem that if a = i?A,o + $a is a well-defined 
self-adjoint operator, defined by its quadratic form [1311221 • Throughout the 
paper unbounded operators will appear only as relatively bounded pertur- 
bations of positive operators in the quadratic form sense, so, in order not 
to complicate arguments and keep the notation simple, we will typically not 
distinguish between operators and corresponding quadratic forms. 

We will assume now for simplicity that A is a cubic volume with periodic 
boundary conditions. Clearly, in this case f^A.o is a non-degenerate ground 
state of i?A,o with a spectral gap: 

#a,o^a,o = 0, -ffA,o|w A e^ A ,o ^ |Ao|l- 

The following result is a perturbation theory for the ground state in the case 
of small a and (3. 

Theorem 1. There exist positive a and (3, depending only on the dimension 
v and the interaction range A , such that if condition holds with these 
a,j3, then: 

1) H\ has a non-degenerate gapped ground state Qa '■ 

Ha^a = Ea^a, 
and for some independent of A positive 7 

H A \ n A en A 

2) There exists a thermodynamic weak* -limit of the ground states VLa '■ 

(AQ A , fi A ) ^> u(A), A e UiakocWa), 
where £>(7Ya) is the algebra of bounded operators in T~Ca- 
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3) There is an exponential decay of correlations in the infinite volume 
ground state uj : for some positive c and e < 1 

WAiAz) -u(A 1 )u(A 2 )\ < c^l+l^le^^^llAillll^H, A, G B{H K ). 

4) If within the allowed range of perturbations the terms <p x ( or the resol- 
vents {h x + 4> x — z)~ l in the case of unbounded perturbations) depend analyt- 
ically on some parameters, then the ground state uj is also weakly* analytic 
in these parameters (i.e. for any local observable A its expectation uj(A) is 
analytic). 

Example 1 (anharmonic quantum crystal model). Let 7i x = L 2 (R d ,dq) 
and 

H = ^(-A x + V 1 (q x )) + \ Yl V ^%)- 

x \x-y\=l 

Suppose that Vi(q) — > +oo as q — > oo. In this case —A + V\ has a discrete 
spectrum with a non-degenerate ground state. Since —A > 0, we see that if 
for some c\ , c 2 

\V 2 (q x ,q y )\ < Ci{Vi{q x ) +V x {q y )) + c 2 , Mq xi q yi 

then for sufficiently small coupling constant A the operator is well-defined 
by the KLMN theorem, and Theorem 1 applies. 

The next theorem extends the perturbation theory to all a G (0, 1) at 
the cost of a slightly more stringent assumption about the perturbation. We 
replace (j2J) with the following stronger assumption: 

^ = 0i r) + ^ b) , (3) 
where <jy- r ' is the "purely relatively bounded" part of the perturbation: 

\<Pt\v,v)\<a\\h x l 2 v\\\ (4) 
and 0^ is the bounded part: 

\\<t>i h) \\<P- (5) 

In particular, (jlj) and (JTJ) imply that (fix Qao+ Xi o = if (f) x ^ is viewed as an 
operator (more precisely, 4> x \v, Q\ 0+x $) = for all v G Dom(fti/ 2 )). 

Theorem 2. For any x > 1 there exists 5 = 5(>c, u, A ) > such that: 
for any a G (0,1), if conditions (GJ)-(Ej) are satisfied with this a and (3 = 
5(1 — a)*( u+1 \ then all conclusions of Theorem 1 hold. 
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Remark. The assumption (jlj) can be somewhat relaxed. In fact, what is 
actually used in the proof of Theorem 2 is not (j3J) but the weaker condition: 
for any I C A 

\Y,^\v,v)\<a\\H]! 2 v\\\ (6) 

This is the condition which we will use when we consider the AKLT model. 
Example 2. Consider a Hamiltonian 

.r 

where A x is a self-adjoint operator on Ha 0+x such that 

A x fl Ao+xfi = 0, A x \ HAo+xe n Ao+xfi > 1. 

Clearly, ® X Q X is a ground state of H with a gap > |A |. We expect that a 
perturbation theory in the sense of Theorems 1,2 holds at least for general 
weak bounded perturbations of H . Theorem 2 shows that this is indeed so 
at least if \\A\\ < oo (A is the operator whose translates A x s are). Indeed, 
consider a finite range perturbation ^2 x ip x with small H^xll- By some rear- 
rangement of terms in H, we may assume without loss of generality that ip x 
acts on H\ 0+x . Now, let A x = h x + 4> x \ where 

h x = mil-fwAo+arefiAo+^o' ^ — A x — \\A\\Pn Ao+x en Ao+Xi0 - 

Here and in the sequel Px stands for the projector onto X. It follows that 
Y2 X h x is a classical Hamiltonian satisfying our assumptions and, by the spec- 
tral gap condition on A, Y^ x </>x is its relatively bounded perturbation so that 
(jU) holds with a = (\\A\\ — 1)/||A|| < 1. We consider now ijj x as (J) x b \ and 
then Theorem 2 applies. 

Now we describe the application of Theorem 2 to the AKLT model. This 
model was introduced by Affleck et al. 0] as the first rigorous example of a 
system in the Haldane phase f |12H 13j. see pQ for a review of the Haldane con- 
jecture). It is a spin-1 chain with the translation-invariant nearest-neighbor 
isotropic interaction 

H = P (2 \Sk + S fc+1 ) = ]T(S fc • S k+1 /2 + (S fc ■ S fe+1 ) 2 /6 + 1/3), 

fcez fcez 

where is a spin-1 vector at site k, and P^ 2 ^(Sfc + S^+i) is the projector 
onto the subspace where + S^+i has total spin 2. The AKLT model has 
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a unique gapped qround state u minimizing the energy of each term in the 
interaction: 

oo(P ( - 2 \S k + S k+1 )) = 

(a frustration-free ground state). The state ui can be described as a valence- 
bond-solid state [S| or a finitely correlated state fOlE]- On a finite chain A 
with periodic boundary conditions the AKLT Hamiltonian H\ has a unique 
frustration-free ground state. 

Let $ = Ylk be any translation-invariant finite range interaction on 
the spin-1 chain. We consider the perturbed AKLT model H + <3>, starting, 
as before, with periodic finite chains A. We prove 

Theorem 3. If \\<fik\\ < ft, with some (3 depending on the range of $ ; then 
all conclusions of Theorem 1 hold for the perturbed AKLT model H + $. 

The main point of the proof is that at large scale the AKLT model is a 
relatively bounded perturbation of a classical model. This enables us to use 
Theorem 2. Though in this paper we restrict our attention to the AKLT 
model only, this property is definitely more general; one can expect some 
form of it to be generic to non-critical gapped spin systems. 

2 Proof of Theorem 1 

We follow the standard approach and approximate the ground state with low- 
temperature states. After time discretization we obtain a cluster expansion, 
which identifies the model with a low density hard-core gas of excited regions 
on the space-time lattice J7j. After that all conclusions of Theorem 1 follow 
in a usual way. Our exposition is, however, rather different technically: we 
derive necessary cluster estimates using the Schwarz lemma and resolvent 
expansions instead of the Feynman-Kac formula. 

We begin by proving that H A has a gapped ground state. Fix some to > 
and consider the expectation 

z NA = ((e- toH ^ N n Afi ,n A>0 ), (7) 

at large AT G N. If Q\ is a non-degenerate ground state of if a with the energy 
E A and a spectral gap > 7, then 

and hence, if Oa >0 ) 7^ 0, 

In Z N>A = 2 In I (ft A , fi A , ) | - t E A N + 0(e^ N ). 
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Conversely, if we show that for some constants ai, a 2 , 

\nZ NA = ai + a 2 N + 0(e- a3N ), (8) 

with 03 > 0, this will imply that in the cyclic subspace generated by Q^,o 
the operator Ha has a gapped ground state. We will argue later that the 
asymptotic (|SJ) holds, with the same a 2 and 03, if we add a small perturbation 
to f^A.o in 0, so Ha has a gapped ground state in the whole space TLa- The 
non-degeneracy of the ground state can be deduced by a continuity argument 
from the non- degeneracy of the ground state in the non-perturbed system. 
We begin proving (JHJ) by writing the identity 



-t H A 



IcA 

where 

T AI = ^(_i)l J l-l- / l e -*o(-H'A,o+E :cSJ ^). 
Jci 

Here the operator Ha,j = Ha$ + ^2 x &j ^x * s defined by the KLMN theo- 
rem, like if a- When formally Trotter or Duhamel expanded, Taj is, by an 
inclusion-exclusion argument, the contribution to the total evolution from 
the perturbation of the classical evolution containing terms <fi x with x G / 
(see (IZj). We do not explicitly use these expansions, however. Since all 
non-commutative terms <p x in Taj lie in A/ = U z6 /(Ao + x), we can write 

T AJ = T; e -* oHA \ A /<°, (9) 

where H A \ Al ,o = ExeA:(A o + a; )nA / =0 and T 'i is defined as T A j with H Afi 
replaced by X^:(Ao+x)nA 7 ^0 h x = #A,o - #a\a,,o- For any Ai C A we will 
denote its neighborhood U a; .(A o+x ) n A 1 ^0(A o + x) by K x , so that T[ acts on 

n A, 

Lemma 1. \\T{\\ < (2ae to(3/a ) w . 

Proof. For some J C I, let zj = (z xi , . . . , z x .j.), Xk G J, be a complex vector 
and consider the operator-valued function 

Hj(zj) = ^ h x + ^2 z x<f>x- 

x:(A o +x)nA I ^0 x£j 

If all \z x \ < 1/a, then, by 0, the quadratic form ^2xeJ Zx< ^ x * s bounded 
relative to Ylx-(A o +x)nA I ^0 hx, with a relative bound < 1 : 



^2z x (f) x (v,v)\ < max \z x \ a ( ^ K 
xeJ x -.(A o + x )r\A I ^0 



> 1/2 


2 


) V 


+ max \z x \\ 



(10) 



7 



Therefore Hj(zj) is an analytic family of m-sectorial operators on 7i A for 
z j ^ {\ z x\ < l/a\x G J} (see [Ej). Since by (fTUj) the numerical range 
{(Hj(zj)v,v)\v G Dom (Hj(zj)), \\v\\ = 1} of these operators lies in the 
half-plane {Rez > \J\(3/a}, it follows from the Hille-Yosida theorem that 

^ e -t Hj(zj) || < e t Q \.J\f3/a^ ^-Q 

Now we consider the operator-valued function 

Tl ( Zl ) = J2(-l) lIHJl e- toHj(zj \ 
Jci 

where zj is a restriction of zj to J. The function Ti(zi) is analytic in {\z x \ < 
l/a\x G J} and, by QJ, ||T/(z/)|| < 2lV°l / l /3 / a . Note that if z x = for 
some x G /, then Tj(zi) = because in this case the terms J \ {x} and 
J U {x} make opposite contribution. Finally, T[ appearing in Q is the value 
of Tj(zi) at zj = (1,1, ... , 1). Now we use a many-dimensional version of the 
Schwarz lemma. 

Lemma 2. Let f(zi) be an analytic function in {\z x \ < a\x G /} and 
\f(zi)\ < M for all zj. Suppose that if z x = for some x, then f(zj) = 0. 
Then \f( Zl )\ <Ma-W U xe i\^\- 

This lemma follows by induction from the usual one- dimensional Schwarz 
lemma. Applying it to Ti(zi), we obtain the desired estimate. □ 

By expanding e~ toHA in Ta,j we have isolated the regions with non- 
classical evolution; to obtain the final cluster expansion we need to iso- 
late in addition regions with classically evolving excited states. Denote 
A j = U ze /(A + x), and also 



H' x = H x Q x , H' A = ® x< z Al Tt' x 



and write T A i as 



T A j - T A>I ^2 p H'jPn (A \ Al )\j,o- 

JCA\A/ 

Now define a configuration C as a sequence {(Ik, Jk)\k = !>•••> N}, where 
Jit C A \ Aj h ; it follows that 



'N,A 



where 



w(c) = (j[ (t^Ph'.p^^) n A , , O (12) 



with the time-ordered product YIk=i Ak = An ■ ■ ■ A\ 
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Lemma 3. \w(C)\ < Uk=i [i^ae^'^^e 



fe/3/a^l4L-<o(|Jfc|-|Ao| 3 |/ fc |) 



Proof. We estimate the norm of the operator in round brackets in (j!2|) . By 
Lemma 1, ||Tf|| < (2ae' o/3 / Q )l 7 l. Next, if J C. A \ A/, then | J n (A \ Aj)| > 
| J| — |A | 2 |A/| > | J\ — |A | 3 |/|. Any x G A\A/ belongs to |A | sets of the form 
A + y, these sets don't overlap with Aj and all contain |A | sites; therefore 
by the spectral gap assumption about h x 

H A \ Alfi > |Jn(A\Aj)|l> (|J|-|A | 3 |/|)1. 

It follows that the norm of the expression in brackets in (fT2"J) does not exceed 
(2ae* o/3/a ) |/fcl e"* o(|Jfel " |Ao|3|/fcl) , which implies the desired estimate. □ 

Now for a configuration C we define its support supp C C {0,1,..., N} x 
A as the set 

{(k, x)\k = 0, . . . , x e A h U A h+1 U J k U J k+1 } (13) 

(with A/ = A/ JV+1 = J = Jn+i = We say that configurations are disjoint 
if they have disjoint supports. If C\ and C 2 are disjoint, we naturally define 
their union C = C\ U C 2 as the configuration with I k = U Ij® , J k = 

Lemma 4. If C± and C 2 are disjoint, then w[C\ U C%) = w(Ci)w(C2)- 

Proof. For the configuration C = C\ U C 2 and any n — 1, . . . , N consider the 
vector 

n 



k=l 

so that w(C) = (vn,&a,o)- We have v n = u n ® ^a\(Aj uj„),o with some 
u n e Haj uj n - Analogously, we can define Vn , Vn , U n , U n for C U C 2 . Let 
K n = Aj n U J n and similarly define 1Q , K n 2 ^ for Ci, C2. Since supp C\ and 
supp C2 are disjoint, it follows in particular that K n 1 ^ and K n 2 ' are disjoint, so 
that H.K n = 7~t„(i) ®7i „(2) . We will prove by induction that u n = Un^<S>u n , at 

^A n i\ n 

n = N this implies the desired equality w(C\ U C 2 ) = w(Ci)w(C 2 ). Suppose 
that u n -\ = Un-i ® u n-i- Note that we have in T~Ck UK _ the equality 

U n ® n^^^,, = T^Pn^^^^ ^Un-l <g> fi(^ ttU jf B _ 1 )\j r - n _ li o)> ( 14 ) 

where T' Kn is defined as T A ,/ with H A , replaced by J2 X :(A o + x )nK n ^0 h x = 
Ha,o - HA\K n ,o- (ffU) holds because e~ to( - HA -°~ HA \ K ^ acts trivially on the 
ground state. By the disjointness, the objects in (j!4)l factor into products of 
respective objects for Ci,C 2 , which proves the inductive step. □ 
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A polymer \ is a connected configuration (i.e., which is not a union of 
two configurations with disjoint supports). We have 

EJV.A " 
Q1SJ. XlvXn J - J - 

k=l 

where summation is over all disjoint collections of polymers in {1, ... , N} x A. 
This is the desired polymer expansion. By Lemma 3, for any e > we can 
choose to large and then a, (3 small so that w(x) < gl su PPxl. \ standard 
combinatorial argument shows that the number of polymers with | supp x\ — 
n containing a given point does not exceed c n for some c = c(u, A ). Now all 
conclusions of Theorem 1 follow from standard results on cluster expansions 
1201 EH E3 E3 Ej, and we will be very sketchy. We define a cluster X 
as a connected collection of polymers xi, ■ ■ ■ > Xk with positive multiplicities 
n>i, . . . , nk- Let G(X) be a graph with rii + . . .+nk vertices, corresponding to 
these polymers, and a line between two vertices drawn if the corresponding 
polymers intersect. Let Gi < G(X) stand for a connected subgraph Gi 
containing all vertices of G(X), and l(Gi) be the number of lines in G\. 
Then the weight of the cluster X is defined as 

w(X) = (n 1 \---n k \r 1 w(XiT 1 ---w( X kT k ]T (-l) 1 ^ . 

Gi<G(X) 

It follows that 

A — ^disj. Xi>--Xn - LJ - A — 'X 
k=l 

with the absolutely convergent series on the r.h.s. (see jHl EH] for recent 
simple proofs). Let l(X) be the time length of a cluster; shifting clusters in 
time, we write 

w(X) = Y (N — l(X))w(X) 

X V ' ^X:l(X)<N K V " V ' 

= -E JC KX)w(X) + N J2 x w(X) + Y X:l{x)>N (m-N)w(X), 

where ^ t_0 ' A i§ the sum over clusters starting at t — 0. By the cluster 
estimate, all series in the r.h.s. converge absolutely, and the last term is 
0(e N ) because summation is over clusters with length > N. Comparing this 
with (jHJ), we identify ai as — Y^x °' A KX)w(X), a 2 as Y^x °' A W {X), and a 3 as 
— lne. This e does not depend on A, so the spectral gap estimate is volume- 
independent. To complete the proof of 1), we consider the changes in the 
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asymptotic (jHJ) when tt A) o is replaced by £l\,o + v with small v. This replace- 
ment adds new polymers Xv, arising from the new terms ((e~ toHA ) N v,Q^ ), 
((e~ toHA ) N n A ,o,v) and ((e" toffA )\i;). The support of Xv contains {0} x A 
or {N} x A, or both. For v small enough the estimate |u>(x)| < el su PP*l re- 
mains valid for Xv The expansion for lnZjv,A is modified by adding clusters 
containing the new polymers Xv Such clusters touch the boundary of the 
time segment {0, . . . , N}, and hence their contribution is c + 0(e ). This 
completes the proof of 1). To show that 2) holds, one writes 

(A{l A ,tt A )= lim Z^ A (Ae- toNHA n Afi ,e- toNHA n Afi ), 

using the fact that (Q At o,fl A ) ^ 0. If A acts on H Aa , one introduces new 
polymers xa with the support containing {0} x and the weight calculated 
using A inserted in the Oth layer; one has < el su PPXA|-|A A |||^||_ it 

follows that (AQ A , Cl A ) = Y^x w (X), where the sum is over clusters in Z x A, 
containing one polymer xa with multiplicity 1. As A /* Z", this expression 
tends to the absolutely convergent sum over polymers inZxZ", which proves 
2). 3) follows from the fact that the truncated correlation on the l.h.s. equals 
^2 X W (X) over clusters containing either a polymer with the support con- 
taining {0} x A x and {0} x A 2 , or two polymers XAi, Xa 2 - Finally, 4) follows 
because if <p varies analytically, then the cluster expansion does too and is 
convergent as long as the estimate (J2J) holds. 



3 Proof of Theorem 2 

Following [TZj, in order to extend the perturbation theory to a close to 1 
we use a scaling transformation. We group lattice sites in cubic blocks b x 
of linear size I, so that the initial cubic volume A is transformed into cubic 
volume A whose sites x are these blocks (we assume that the size of A is 
a multiple of I, but one can consider general cubic volumes too by taking 
blocks of different sizes). For any x G A, let 

T~t x = ®y£b x 7~(-yi = ®y£b x ^yi 7~i x = © 

and also 

for I C A. Suppose that I > diam(A ). We can then view the interaction 
$a as the sum J^xgX^xj here 4> x ac ^ s on ^a +z> where A = {0, 1} U , and is 
defined as 

2/:Ao-H/eU z£ T +x bz 
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Here c y = \{x : A + y G U 2g ^ o+x 6 z }| is the number of (p x where (f> y appears. 

We define analogously for (j)x\(f^x' appearing in the decomposition 

0. Next for any / c A we define 

x£l 

— i r \ — — (r) — (b) 

and similarly for $},<£>}, so that H A = H Afi + $j + $ j . Following the 
proof of Theorem 1, we write 

e~ toHA = J2_ T ^ 

I,J,KcA 

where 

r w = £^(-i)i'H'.i(-^H J .i e -«^<«S5, 

IlCJ JiCJ 

We call a sequence C = {(Ik, Jk, Kk), k = 1,...,N} a configuration and 
assign to it the weight w(C) = (H k=1 Ti k ,j k ,K k tt& )0 , ^a,o)- Let ^l,---,^ 
stand for the coordinates of the site x G A; for any set I C A we define its 
neighborhood 

/ = {x\3y G / such that \x k — yu\ < 1, k = 1, . . . , z/}. 

Let A/ = U xg /(A + x). Similarly to (fT3|) . we define a configuration's support 

as 

{(k, x)\k = 0, . . . , A^; x G % h U 5 /fc+1 U X Jfe U A Jfc+1 U K k U AT fc+ i}- 

An analog of Lemma 4 on factorization of weights is immediate. Therefore 
the only thing that needs to be proved is an exponential bound for the weight: 

\w{C)\ < e l suppC l 

with sufficiently small e; after that the conclusion of the theorem follows like 
in the previous section. Clearly, in order to have this bound it suffices to 
show that for any e one can choose to, I and (3 so that 

wt^kW < c i'i+u+in (i6) 

The remainder of this section is a proof of this claim. Specifically, we will 
show that for a close to 1 one can achieve this by choosing 

*o = (l-a)-*, 1= f(l -<*)-"], P = 5(l-a)^ +1 \ (17) 
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where ["•] is the integer part; here x is any fixed constant > 1, and 6 = 
S(k, v, A ) a sufficiently small constant. The strategy of the proof is as fol- 
lows. We will obtain three different bounds for HT^^H, suitable when the 
contribution to I U J U K from either J, J or K is large enough. 

Case 1. The first bound relies on the smallness of the bounded part (j) 6 
of the perturbation, and is used when J is large. In this case we use again 

the Schwarz lemma. In the definition of Ti,j,k we replace YlxeJ *f*x with 

Xlx-eJi Zx( ^x ■> where z x G C, \z x \ < a, with some a > 1. Let a = (£o||0^||) - \ 

then by definition of ^ and from (|17|) 

a > (to(2Z)1l0 (6) ||)- 1 > (to(2/)^)- 1 > 2- v S~ 1 > 1 
if 5 < 2~ u . Using the Schwarz lemma with this a, we find that 

\\Tjj, k \\ < 2 |/| (2et p (6) ||)l J l < 2l / l(2^ +1 e5)l J l (18) 

Case 2. The second bound relies on the contractiveness of the classical 
evolution in excited regions and is used when K is large. We will estimate 

the norm of e 01 A, ° 7 i 3 ^'Pwj> ^ • We begin by writing 

Kk®&a\k,o to j to 

e -t (H A , 0+ 3> h +$Ji ) = J e -t z Rzd ^ (19) 

where R z is the resolvent (.Ha,o + + ^j? — z ) ^ an d T is a contour in 

the complex plane going around the spectrum of H^ + $^ + $^ ; we will 
specify T below. We will use the expansion of the resolvent 

oo 

R Z = Q Z (J2 F *)Q*> ( 2 °) 

k=0 

where 

F z = -(H A , +^ - z)- 1 ^(H A>0 + ^ ) 1 - zY 1 ' 2 

and 

The operators Q z are well-defined and bounded for z in the resolvent set of 

H\ + 3?^ j i-e., at least in C \ [— ||<&^ ||, +oo). We now estimate the norm 
of F K : 
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|-"A,0 + - 2! | |-"A,0 + - 2; I ' I 



sup 

«eft x \{0} 



lull 2 



sup 



»6Dom(|H A>0 +i™-z|V2) ||l#A,0 + - z| 1/2 t>|| 2 

< SUp ^ 

J ,6Dom(|if A ,o+*Si-z| 1/2 ) IH^A.O + ^J! ~ ^ 1/2 t>|| 2 



< sup 

Ae Spec(ff A ,0+*£ ) ) 



«(A+||*S?| 



A 



We will be interested in those z where \\F Z \\ < y/a. By the above bound, a 
sufficient condition for that is 

z£ |J [zec\\z-X\<^ + (21) 

A6Spec(flA,o+*^) 

Since Spec (i?A,o + C [— ||$j ||, +00), the above union lies in the sector 

{z E C : I arg(z + < arcsiny^}, (22) 

so if we choose z outside this sector we have H-F^H < \/a and in particular 

a k ' 2 

\\Q z F k M < z^-. (23) 

dist (z, Spec(i/A,o + $j;)) 

We will now show that, furthermore, one can enlarge the domain of z 
where a bound of the above type holds, if one applies the operator on the 
l.h.s. to vectors from Ti K <g> fij^p with K large compared to J. Precisely, 
let 

n = \{\K\- 6"| J\)/T\ 

This n is a lower bound for the maximal number of sites in a subset K\ C 
K such that the neighborhoods {x} of points x G K\ are separated from 
each other and from Aj by at least two A-lattice spacings: choose the first 
such x outside the 2- neighborhood of Aj, then the second outside the 2- 
neighborhood of Aj unioned with the 3-neighborhood of the first x, etc. We 
assume that n > 0. Next, let 

m — \lj diam (Aq)] , m > 0. 
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For any a define lA a as the union of circles standing in (|21|). but with A 
running over [a, +00). Now, suppose that \k/m\ = r with r < n; we claim 
then that 

11^*^- Jl < — — ^ „^).. (24) 



dist( 2 ,[|A |(n-r)- ||$^||, +00)) 



*■ |A |(n-r)-||*^|| 

Indeed, let Q[ a ,+oo) stand for the spectral subspace of H^ q + $^ , cor- 
responding to the interval [a, +00). Let K x be chosen as above, \K X \ = n. 

Then the subspace Tt' Kl ®T~Cj\Ki is an invariant subspace of H\ t0 + &^ and it 
lies in Q n , . ..^m,. The sites 2 £ .fTi (and therefore their neighborhoods 

- / [|A |rHI$j 1 'j|,+oo) V to 

{x} too) are in excited states. When we apply powers of F z to vectors from 
TC Kl <8> ^a\a"i' we nee d a ^ least m = [// diam (A )] powers to remove the 

excitation from a neighborhood {x} of a given a; £ Ki, because the block b x 
has to be connected to A\U y£ ^b y by supports of the elementary interactions 

4>w\w £ A, of which $^ is composed. Hence to remove excitations from r 
neighborhoods we need at least mr powers of F z . Therefore if k < mr with 
r < n, then 

F Z H K ® Qa\^,0 C ^[|Ao|(n-r)-||l^||,+oo) 



(here we use 7~tK®^A\K,o c ®^A\ifi)- ^ follows that in the case at hand 

— ( b ) — ( r ) 

we can replace H^ + $ j and the quadratic form $ 7 with their restrictions 

to ^n Ao |/ n _ r N_,|^C(')|| +0O y the argument leading to then yields (J23)- 

Now we can specify the integration contour T. It will depend on the term 
k in the expansion ()2Uj) through r = [fc/m] . Let s = n — r. For s = 0, . . . , n, 
let 



z s ,± = (1 - V5)|Ao|a - ± z-^|A | S - ^ 1 

Vl — a 

and let the contour r s consist of the segment [z S) -, z s ,+\ and the two rays 
{z\ arg(z — z S) ±) = ± arcsin y/a}, so that it is the boundary of the truncated 
sector (J22J) (of the sector itself in the case s = 0) shifted by t$ 1 to the left. 
This contour is just a convenient for calculations, piecewise-linear approxi- 
mation of the boundary of U. K . ,,^(6),,. We choose the contour in this way 

because we need it to lie as far to the right as possible due to the factor e~ toZ 
in (fTTjj) . but still outside of ^| Ao | g _||$( 1, )||j so that the resolvent estimate can be 

used. By slightly shifting it to the left we avoid the possible singularity in 
the denominator in (|24ft. 
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|Ao|s-||*Sil 



We have T S C\U,, . ,,■=&),. = 0, so we can write 

n K<® ll A\K,0 

f (r+l)m— 1 



r=0 Jr "- r ^ fc=rm ' 



v fc=nm 

We estimate now this expression using | J e~ toz f(z)dz\ < J e _ioRe2 |/(2;)||(iz| 

and the bound fl2IJ). We have dist (T s , [|Ao|s-||$$J||, +oo)) > y/a. It follows 
that 



< 



e *oll4ill+i 



n-l 



2 

V^V \/ato 



< 



e <o||*Sill+l 



(n + 1) (|A |rj H — ■=— ) (max{a 2 



-to(l-v^)|Ao|\\ n 



})». 



27TV1 - \/ato 
Some calculation now shows that iito,l,P are defined as in (|T7|) with x > 1, 
then for any e > we have 

lle-M^o+^+^P-w - II < e *oll ¥ Jill e « < e 2 " sw e n 

when a > a for some a = a (e, x, u, A ) < 1. It follows that 

\\Ti,j,k\\ < 2^+IV^ V. (25) 
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Case 3. The third bound is used when / is large. We will bound 
|| ^ 7iC7 (-l)l^l-l- f il e - t °™+*7 1 + * J i ) R# K9TiAKo \\. Like in case 2, we represent 

where R z ,h = (-Ha.o + + ^j? — z )~ l an d T is the shifted boundary of the 
sector defined in case 2. We again expand R z ,h = Qz(J2T=o^zh)Qz- Now, 
let 

n'=\(\I\-V\J\-V\K\)/Vl 

This n' is a lower bound for the maximal number of sites in a subset J 3 cJ 

such that the neighborhoods {A + x} of points x G I2 are separated from 
each other, from Aj and from K. Let m = [7/ diam (Ao)] as before. We 
claim that 

E(-i) |i| -""^ 1 o^h^, = » 

hci 

— (r) 

if fc < mn. Indeed, if in the above sum we expand F Zj / 1 = —Q z ^2 xeIl ^ 
then by the inclusion-exclusion principle it becomes 

k 

(-i) fe E(n(^^))^ p ^% VK0 > 

where summation is over sequences X\, . . . ,Xk of points from J, containing 
each point of I at least once. In particular, each sequence contains all points 
of I2 defined above. Each term in this sum is 0. Indeed, if we further expand 

each (jx^ in (jffi as in (fT5|) and consider the set Up =1 (A + y p ) for each of the 
resulting sequences, then, if k < mn', this set will have at least one connected 

component contained in the set {A + x} for some x 6 I^- By the choice of 
I2, the operator (f>ty from this component which acts first acts on the ground 
state. Since 4>^Q\ 0+y = 0, this implies our claim. 
It follows that 

\\Y(-i)W-m e ^ s ^^^p# ^ 

\\^ K ' W K ®%K 
hCl 

/oo 

jiu -° k=mn' 

7rto\/«(l — a)(l — y/a) 
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Again, for any e > 0, if a is sufficiently close to 1, then this expression does 
not exceed 2' / 'e 2 ' / ' s ' Jl 'e n ', which implies 




< 2 \I\+\J\ e ^S\J\ e n' 



(26) 



in case 3. 



The desired bound ()16|) now follows from bounds obtained in the three 
cases. If / is sufficiently large compared to J and K, then one uses the bound 
(|2T)J) . Otherwise, and if K is sufficiently large compared to J, one uses (|23|) . 
In the remaining case one uses (fTSj) . 

4 Proof of Theorem 3 

We summarize some known facts about the AKLT model which we will need. 
Denote by H\ and the AKLT Hamiltonians on a finite chain A with 
periodic and free boundary conditions, respectively. The Hamiltonian has 
a four-dimensional subspace Qa of frustration-free ground states. Using the 
valence-bond-solid representation, one can choose a (non-orthogonal) basis 
^A-,ab, a,b = 1, 2, in with the following properties: 
1) For two adjacent finite chains A 1; A 2 



Then g A = 1 + 0(3"I A I) as |A| -> oo. 

We will also use the fact that the operators have a uniformly bounded 
away from spectral gap, i.e. for some 7 > we have > 7(1 — G\) for 
all A, where Ga is the projector onto Q\. 

All the above facts were proved in see also |51 EH ITT1 IT51 IT?! 1131 124] for 
various refinements. 

We show now that at large scale the AKLT model is a perturbation of 
a non-interacting model. Like in the previous section, we group the sites of 
the cyclic chain A in blocks Ai, A2, . . . , A n with |A^| = I, k = 1, . . . , n. We 
will specify I later. We write 



^AiUA 2 ;afc — ^Ai;al ® ^A 2 ;26 — ^Ai;a2 ® ^A 2 ;lb- 



The unique ground state of H\ is given by £3a ; i2 — ^a ; 2i- 
2) Let g\ be the 4x4 Gram matrix of the basis fl\ ;a b'- 



(9 a) ab,cd (^A;ab, ^A;cd)- 
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k=l 
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where 

H^k+i = H{J2 + P®(S kl + S fem ) + H{ k+ J2 

= H{ kUAk j2 + pW(S kl + S kl+1 )/2 

(with the convention n+ 1 = 1). Clearly, Ker (H k)k+ i) = Ker (H AkUAk+ ^) = 
^A fc uA fe+1 and 

H k , k+ i > H{ kUAk j2 > 7/2(1 - G AfcUAfe+1 ). (27) 

Now, an important role is played by the asymptotic commutativity of the 
projectors G\ k u\ k+1 , which we utilize as follows. For each k we orthogonalize 
the basis ^A k - a b- 

2 

^A k ;ab = (9Al /2 )ab,cd^A k ;cd 
c,d=l 

and next define 

^A fe UA fe+1 ;afe = ^A k ;al ® ^A k+1 ;2b ~ ^A k ;a2 ® ^A fe+1 ;l&- 

Denote by G" k fe+1 the projector onto the four-dimensional subspace spanned 
by fiJL uA fc+i; ab i n ^AfcUA fe+1 - A straightforward calculation shows then that 
G" k _ x k commutes with G kk+V At the same time, by the property 2) above 

\\ G k,k+i ~ G A k uA k+1 1| = 0(e l ) (28) 

with some e < 1. Now we use the following abstract observation. 

Lemma 5. Let Hi,H 2 ,H^ be three finite- dimensional Hilbert spaces, and 
Hi and H 2 be two commuting self-adjoint operators acting on Hi ® 7^3 and 
7^3 ®H 2 , respectively. Then there exists a decomposition 

such that for each s Hi ® H\ s ® H\ s and H\ s ® 7Y| s <S> H 2 are invari- 
ant subspaces of Hi and H 2 , respectively, and, furthermore, the restriction 
Hi\ ni(gjn i ^ H 2 is an operator acting only on Hi ®H\ S , and the restriction 
H 2 \ H i ^2 0W2 is an operator acting only on H\ s ® H 2 . 

Proof. Decompose Hi = £\ Hu i ® Hi 3 i , where Hn^ are linearly indepen- 
dent operators on Hi, and Hi 3i are operators on H 3 . Consider the algebra 
Ai C B{H-$) generated by the operators H 13t i and the unity. Ai is a von 
Neumann algebra (closed under taking adjoints). Similarly, we decompose 
H 2 = J2i^23,i <8> H 2 2,i, where H 22ji e B(H 2 ) are linearly independent and 
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^23,% £ B(7i 3 ). By the commutativity assumption the operators H 2 ^^ lie in 
the commutant A[ of A\. But it follows from the well-known classification 
of finite-dimensional von Neumann algebras (see e.g. [Zj) that there exists a 
decomposition 

n, AHl : : Wg,.) (29) 

such that 

A = ® s (B(Hl,) l wg J, ^ = © s (l wis B(Wl,J). 

Therefore (|29jl is the desired decomposition. □ 

We apply this observation to the Hilbert spaces 'H\ k _ l ,'Hh k ,'H\ k+1 and 
operators G k _ y k , G' k fe+1 . The relevant decomposition then is 

H Ak =Ft®F!®(H Ah egA h ). (30) 



Here T\,T k are two-dimensional Hilbert spaces such that T\® = (?A fc . 
One can choose orthonormal bases {v\ a } a= i^ and {^ fe }&=i,2 in and T k 
so that ^A fc ;afe = v l,a® v l,b- The subspace 7^A fe _ 1 £?A fe is invariant for G k _ lk , 
and in this subspace G" k _ lk essentially acts only on T\; similarly, 

Qk k ®^A fc+ i is invariant for G" k k+1 , and it acts there only on ^ 7 |®HA fc+1 - The 
subspaces H.A k _ 1 (7YA fc £7aJ and (Ti.A k C/aJ %v fc+1 he in the kernels 
of fc , fe+1 , respectively (hence we don't factor H.A k QA k in P0"| ). 

If we use the decomposition (|H0J) for two neighboring blocks k and k + 1, 
we get a decomposition of HA fe ^A fc+1 as a direct sum of four subspaces: 

HA k ®H Ak+1 = (^A fc 0£A fc )0(H Afc+1 0£A fc+1 ) 

© (^0^)0^+10^+1 
© ^0^0(HA fc+1 0^A fc+1 ) 
© ^0^ 0^+1 0^+i- 

The first three subspaces lie in the kernel of G kk+1 , whereas in the fourth 
G" k k+1 acts as the projector onto the vector 

Vk,k+1 = Vfc,i ffe+1,2 - 4,2 © vl+i,i ( 31 ) 

in the space T k 

In order to get a classical model in the sense of Introduction we introduce 
additional Hilbert spaces T k ,T k , with dimjF| = dim^ = 3'/ 2 — 2 (assuming 
that / is even), so that 

^A fc = (^©^ 3 )0(^©^ 4 ). 
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Now define new Hilbert spaces Tl k)k +i by 

H k , k+1 = {Tl © Tt) ® (fUx © ^fc+i)- 

The initial scaled spin chain with sites indexed by k and Hilbert spaces 7i.A k 
assigned to k is then equivalent to the chain with sites indexed by pairs 
(k,k + 1) and Hilbert spaces TLk,k+i assigned to the new sites (k,k + 1). 
Let h kjk +i be the projector in H k , k +i onto the orthogonal complement to the 
vector v k> k+i introduced in pijl. Consider the operator 

n 

H\ t o = 3/ hk,k+i- 

k=l 

We claim that if I is large enough (independently of n), then this operator 
is the desired classical Hamiltonian, such that the AKLT Hamiltonian H p h is 
its relatively bounded perturbation satisfying assumptions of Theorem 2. 
To prove this, we write as 



:,fc+l' 

k=l k=l 



where 



^k,k+i — C 1 _ G' k \ k+ i)Hk,k+i(l - G'k,k+i) ~ Khk-i,k + h k ,k+i + h k +i,k+2) 
will be the "purely relatively bounded" part of the perturbation, and 

4>k!k+i = H k ,k+i - (1 - Gfe ifc+1 )-fffc,fc+i(l - G'l k+1 ) 
the bounded part. First we estimate ||0fcL_ilh 

II^MH-lll = IK 1 ~ G A k UA k+1 )H k , k+ i(l - G Afe uA fc+1 ) 

_ (! " g aU+i)-^M+i(1 - G"fe,fc+i)ll 
< ||CA fc UA fc+1 - G kk+1 \\(\\H ktk+1 (l - G Ak uA h+1 )\\ 

+ \\(l-Gl k+l )H ktk+1 \\) 
= 0(le l ), (32) 

by ()28j) and because Hi^^iH < /. Now we analyze the term X^fc^fcjfc+i- 

We 

claim that the condition © holds with a = —(1 —7/(6/) +0(e )), uniformly 
for all J C {1,2,..., n}. Indeed, note first that 

h k , k +i < 1 — G k k+l < h k -\ jk + h k)k+ \ + h k+ \ )k+ 2- (33) 
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It follows from the right inequality that <j> k k+l < 0. Therefore the maximum 
over I in (0) is attained when I = A. By (J27|) . (J2%|) and the left inequality in 

(! - G l,k+i) H k,k+i(l - G'k,k+i) 

> 7/2(1 - G'l k+1 )(l - G AfeUAfe+1 )(l - G' k \ k+1 ) 

>( 7 /2 + 0( e '))(l-G^ +1 ) 

> (7/2 + 0{e l ))h k ,k+l 

and hence 

n n 

-l,fe + ^fc,fc+l + ^fc+l,fc+2)) 

fc=l fc=l 

= -(l- 7 /(6/) + 0(e i ))^A,o. 

Since < 0, this proves our claim about relative boundedness with 

a = 1 -7/(6/) + 0{e l ). 

Now we apply Theorem 2. We have (1 - a)^ 1+u) = ( 7 /(6/) + 0{e 1 )) 2 *. 
On the other hand, by the bounded part of the perturbation is 0(le l ), 
which is asymptotically less than (7/(6/) + 0(e l )) 2>c . Therefore for I large 
enough H P K is a relatively bounded perturbation of i?A,o so that Theorem 2 
is applicable. The conclusion of Theorem 3 follows now from Theorem 2, 
because a sufficiently weak perturbation of the AKLT model remains within 
the range of perturbations of H^ Q) where Theorem 2 is applicable. 
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